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Determining the waveguide conductivity in a hyperbolic 
equation from a single measurement on the lateral boundary 


1 Michel Cristofol, 1 2 Shumin Li, 3 Eric Soccorsi 

Abstract 

We consider the multidimensional inverse problem of determining the conductivity 
coefficient of a hyperbolic equation in an infinite cylindrical domain, from a single boundary 
observation of the solution. We prove Holder stability with the aid of a Carleman estimate 
specifically designed for hyperbolic waveguides. 
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1 Statement of the problem and results 

1.1 Introduction 

The present paper deals with the inverse problem of determining the time-independent isotropic 
conductivity coefficient c : P —> R appearing in the hyperbolic partial differential equation 
<9f — V • cV = 0, where P := uj x R is an infinite cylindrical domain whose cross section ui is 
a bounded open subset of M n_1 , n > 2. Namely, t > 0 being arbitrarily fixed, we seek Holder 
stability in the identification of c in := uj x (—£,£) from the observation of u on the lateral 
boundary Yj j := du x (— L, L) over the course of time (0, T), for L > £ and T > 0 sufficiently 
large. 

Several stability results in the inverse problem of determining one or several unknown co¬ 
efficients of a hyperbolic equation from a finite number of measurements of the solution are 
available in the mathematics literature [1, 5, 6, 7, 17, 18, 26, 29]. Their derivation relies on 
Bukhgeim-Klibanov’s method [9], which is by means of a Carleman inequality specifically de¬ 
signed for hyperbolic systems. More precisely, [17, 29] study the determination of the zero-th 
order term p : P —> R appearing in dt — A + p = 0, while [1, 7] deal with the identification of 
the speed c : H I in the hyperbolic equation dt — cA = 0 where A = A(x, D x ) is a second 
order differential operator. The case of a principal matrix term in the divergence form, arising 
from anisotropic media, was treated by Bellassoued, Jellali and Yamamoto in [5], using the 
full data (i.e. the measurements are performed on the whole boundary). Using the FBI trans¬ 
form Bellassoued and Yamamoto claimed logarithmic stability in [6] from arbitrarily small 
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boundary observations. Imanuvilov and Yamamoto derived stability results in [18] by means 
of H -1 Carleman inequality, from data observation on subdomains fulfilling specific geometric 
assumptions. In [26] Klibanov and Yamamoto employed a different approach inspired by [25]. 


Similarly, numerous authors have used the Dirichlet-to-Neumann operator to claim stability 
in the determination of unknown coefficients of a hyperbolic equation. We refer to [4, 19, 28] 
for a non-exhaustive list of such references. 

In all the above mentioned papers, the domain was bounded. Recently, the Bukhgeim- 
Klibanov method was adapted to the framework of infinite quantum cylindrical domains in 
[10, 23, 24], Inverse boundary value problems stated in unbounded waveguides were also 
studied in [11, 12, 13] with the help of the Dirichlet-to-Neumann operator. In all the six 
previous articles, the observation is taken on the infinitely extended lateral boundary of the 
waveguide. The approach developed in this paper is rather different in the sense that we 
aim to retrieve any arbitrary bounded subpart of the non-compactly supported conductivity c 
from one data taken on a compact subset of the lateral boundary. This requires that suitable 
smoothness properties of the solution to (1.1) be preliminarily established in the context of 
the unbounded domain D. 

The paper is organized as follows. Section 2 is devoted to the analysis of the direct problem 
associated with the hyperbolic system under study. In Section 3 we prove a global Carleman 
estimate specifically designed for hyperbolic systems in the cylindrical domain fh Finally 
Section 4 contains the analysis of the inverse problem and the proof of the main result. 

1.2 Settings 

1.2.1 Notations 

Throughout this text we write x = (x ', x n ) £ D for every x' := (xq,..., x n _i) £ oj and x n £ R. 

/ \ V 2 

Further, we denote by \y\ := (X)q=i Vj ) the Euclidian norm of y = (y\,... ,y m ) £ R m , 

m £ N*, and we put S n_1 := {x/ = (xq,..., x n -\) £ R™ -1 , \x'\ = l}. 

For the sake of shortness we write dj for d/dxj, j £ N* +1 := {m £ N*, m < n + 1}. For 
convenience the time variable t is sometimes denoted by x n+ \ so that dt = d/dt = d n+ \. We 
set V := (<9i,...,d n ) T , := (<9i, ..., <9 n -i) T and V X;t = (<9i,..., d n , d t ) T . 

For any open subset D of M m , m £ N*, we denote by H P (D) the p -th order Sobolev space 
on D for every p £ N, where H°(D) stands for L 2 (D). We write (■, -) P) d (resp., j| • \\ P) d) for the 
usual scalar product (resp., norm) in H P (D) and we denote by Hq(D) the closure of Cfi°(D) 
in the topology of H 1 (D) 

Finally, for d > 0 we put := uj x (—d,d), Qd := x (0,T), := du x (—d, d) and 

Srf := duj x (—d,d) x (0, T). 

1.2.2 Statement of the problem 

We examine the following initial boundary value problem (IBVP in short) 



( 1 . 1 ) 
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with initial conditions {Oq , 6\ ) , where c is the unknown conductivity coefficient we aim to 
retrieve. This is by means of the Bukhgeim-Klibanov method imposing that the solution u to 
(1.1) be sufficiently smooth and appropriately bounded. 

Throughout the entire text we shall suppose that c fulfills the ellipticity condition 

c > c m in fi, (1.2) 

for some positive constant c m . Notice that we may assume, and this will be systematically the 
case in the sequel, without restricting the generality of the foregoing, that c m E (0,1). 

Let us now say a few words on the solution to (1.1). In order to exhibit sufficient conditions 
on the initial conditions (#o,#i) (together with the cross section u and the conductivity c) 
ensuring that the solution to (1.1) is within an appropriate functional class we shall make 
precise further, we need to introduce the self-adjoint operator A = A c , associated with c, 
generated in L 2 {Vl) by the closed sesquilinear form 

qA[u] := ||c 1/2 u||oq = / c(x)\u(x)\ 2 dx, u E D(q A ) := Hq(SY). 

J n 

Evidently, A acts on its domain as —V-cV. Since A is positive in L 2 (fl), by (1.2), the operator 
A 1 / 2 is well defined from the spectral theorem, and D(A 1 / 2 ) = D(q A ) = (Q). For the sake 

of definiteness, we set A 0 := I and D(A°) := L 2 (Q), where I denotes the identity operator in 
L 2 (Q), and we put 

A m / 2 v := A^-^iA^v), v € D{A m ' 2 ) := {v E D(A ( ' m ~ 1 ^ 2 ), A 1 / 2 v E D(^4 (m - 1)/2 )}, 
for each m E N*. It turns out that the linear space D(A m / 2 ) endowed with the scalar product 

m 

( v ’ w )d(a™/ 2 ) : = ^2(A : ’ /2 v,A^ 2 w)o,n, 

3=0 

is Hilbertian, and it is established in Proposition 2.6 that 

D(A m ) = {v E H 2m (n)- v, Av ,..., A p ~ 1 v E Hq(Q)}, m E {p - 1/2,p}, p E N*, (1.3) 

provided du is C 2m and c E 1 / E 2m “ 1 ’ 00 (ll). As a matter of fact we know from Corollary 2.7 for 
any natural number m, that the system (1.1) admits a unique solution 

m+1 

uE f] C k ([0,T}-H m+1 - k {n)), (1.4) 

k =0 

provided the boundary du is C m+1 , the conductivity c E lT m,0O (0;M) fulfills (1.2) and 
(6» 0 ,6»i) E D(A( m+1 )/ 2 ) x D(A m / 2 ). Moreover, if ||c|| — CV/ for some constant c A 0, 

then the solution u to (1.1) satisfies the estimate 

m+l 

y, ll' u llc ,fc ([0,T];// m + 1 - fe (n)) < C (ll^ollm+l,n + 11^1 \\m,Cl) , (1-5) 

k =0 

where C > 0 depends only on T, u and c+f. 
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1.2.3 Admissible conductivity coefficients and initial data 

In order to solve the inverse problem associated with (1.1) we seek solutions belonging to 
n| =3 C fc ([0, T]; H 5 ~ k (Q)). Hence we chose m = 4 in (1.4) and impose on c to be in W A,oc {£l] M) 
and satisfy (1.2). In what follows we note cm a positive constant fulfilling 

ll c llvy 4 ’°°(0) < °m- (1-6) 

Since our strategy is based on a Carleman estimate for the hyperbolic system (1.1), it is also 
required that the condition 

a' ■ V x /c > cio in H, (1.7) 

hold for some a 1 = (a \,..., a n _i) € S n_1 and do > 0. Hence, given 0*, a neighborhood of T 
in M n_1 , and c* G W 4,00 (0 * ft fi;R) satisfying 

c* > c m and a' ■ V^c* > do in 0* n H, (1.8) 

we introduce the set A o„ = {o!, d 0 , c*, c m , cm) of admissible conductivity coefficients as 

:= {c G IH 4,00 (H;]R) obeying (1.2) and (1.6) — (1.7); c = c* in 0* nH}. (1.9) 

Notice that the above choice of m = 4 dictates that (0oA) be taken in 0(A 5 / 2 ) x D(A 2 ), 

which is embedded in i7 5 (H) x H\fl) according to (1.3). Furthermore, it is required by the 
analysis of the inverse problem carried out in this article that 9q be in 1F 3,00 (H) and satisfy 

- a' ■ V X '0 O > rj 0 in w*, (1.10) 

for some r/o > 0 and some open subset cj* in , with C 2 boundary, satisfying 

u \ 0* C w* C w. (1.11) 

Thus, for Mo > 0 such that 

i 

||6»o||iv3.°°(q) + ^2 < -An (1-12) 

j=o 

we dehne the set = Q Ut (a 1 , Mo, rjo) of admissible initial conditions (9o,9i) as 

:= {(00A) G (D(A^ /2 ) n 1F 3 ’°°(H)) x D(A 2 ), fulfilling (1.10) and (1.12)} . (1.13) 

Having introduced all these notations we may now state the main result of this paper. 

1.2.4 Main result 

The following result claims Holder stability in the inverse problem of determining c in Qg, 
where t > 0 is arbitrary, from the knowledge of one boundary measurement of the solution to 
(1.1), performed on Ex for L > t sufficiently large. The corresponding observation is viewed 
as a vector of the Hilbert space 

A(£ l ) := H 4 ( 0, T; L 2 {T l )) D H 3 ( 0, T; H^Tl)), 

endowed with the norm, 

ll' U ll^r(E_ L ) := ll tt llH 4 (0,T;L 2 (r i )) + IMllf3(o ) T;.ff 1 (r I ,)) I U € A(Si). 
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Theorem 1.1. Assume that diu is C 5 and let 0* be a neighborhood of T in 1" 1 . For 
a' = (oi,..., o n _i) € § n_1 , ao > 0, c m £ (0,1), cm > c m and c* 6 IT 4,oo (0* nfl;M) fulfilling 
(1.8), pick Cj, j = 1,2, in Ao* (a', a 0 , c*, Cm, cm), defined by (1.9). Further, given Mo > 0, 
7/0 > 0 and an open subset w* C obeying (1.11), chose do £ @ U] *(a', Mo,r]o), defined in 

(1.13), and 9\ = 0. 

Then for any l > 0 we may find L> l and T > 0, such that the f]° k=0 C k ([0 1 T], H 5 ~ k (Q))- 
solution Uj, j = 1,2, to (1.1) associated with {6o,0{), where Cj is substituted for c, satisfies 

||ci - C2\\ H i(n e ) < C ||ki - ■ 

Here C > 0 and k £ (0,1) are two constants depending only on u, £, Mq, t/o, a!, ao, c*, c m 
and cm- 

We stress out that the measurement of the observation data is performed on Tl and not on 
the whole boundary 8Ql- 

2 Analysis of the direct problem 

In this section we establish existence and uniqueness results as well as regularity properties, 
for the solution to hyperbolic (l.l)-like IBVP systems. The corresponding results are similar 
to the ones obtained for hyperbolic equations in bounded domains (see e. g. [15, Sect. 7.2, 
Theorem 6]) but since P is unbounded here, they cannot be derived from them. 

2.1 Existence and uniqueness result 

With reference to (1.1) we consider the boundary value problem 

f d?v + Av = f in Q , , 

\ v(-,0) = g , d t v(-,0) = h in O, 

where /, g and h are suitable data, and we recall from 4 [27, Sect. 3, Theorem 8.2] the following 
existence and uniqueness result. 

Proposition 2.1. Assume that f £ (7 0 ([0,T]; D(A 0 )), g £ D(A 1//2 ) and h £ D(A°). Then 
there exists a unique solution v to (2.1) such that 

djfv £ 0°([O,T];0(d (1 - fc) / 2 )), fc = 0,l. (2.2) 

Moreover we have the estimate 
i 

sup \\d k v(-,t)\\ D ( A (i-k)/ 2 j < C (j|/||c°([ 0 ,T];.D(A 0 )) + Ibll DiA 1 / 2 ) + INI-D(A 0 )) • (2-3) 

k =o 

4 Upon taki ng H := D(A°) = V := D(A 1/2 ) = H^(Q) and a{t\u,v) := (A 1/2 u, A 1/2 v) 0 ,n = 

f Q c(x)Vu(x) ■ Vv(x)dx for u, v G V and all t £ [0, T\. 
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2.2 Improved regularity 

Proposition 2.2. Assume that f G C' 1 ([0, T]; D(A 0 )), g G 11(A) and h G D(A v t 2 ). Then the 
solution v to (2.1) satisfies 

d?v G C'°([0,T];P(^ 2 - fe V2)) ) k = 0)12) 

and we have the estimate 
2 

E SU P II_D(A( 2_fc )/ 2 ) ^ C ^ll/llc 1 ([0,T];D(A°)) + lll/ll-DfA) + ll^'llo(A 1 / 2 )) • 

fc=0 te[0,T] 

Proof. By differentiating (2.1) with respect to t, we check that w := dtv obeys 

dfw + An; = d t f in Q 

w(-, 0) = h, d t w(-, 0) = /(-, 0) - Ag in Q. 

Since /(-, 0) — Ag is lying in D(A°) then we have d k+1 v = d k w G C°([0,T]; P(A( 1_fe )/ 2 )) for 
k = 0, 1, with the estimate 

i 

E SU P \\^t +lv (^ t )\\D(A( 1 - k )/ 2 ) < C (\\dtf\\c°([ 0 ,T]-,D(A 0 )) + ll^-||_D(A 1 / 2 ) + ||/(‘, 0) - ^IIdCA 0 )) 
fc=0 te[o,T] 

< C (ll/llc ,1 ([0,T];.D(A 0 )) + ||ff||.D(A) + II^IId(AV 2 )) ( 2 -7) 

by Proposition 2.1. Further, as Au = / — d 2 u from the first line in (2.1), we get that v G 
C°([0,T]; D(A)), and that ||'u(-, i) || z?(A) is majorized by the right hand side of (2.7), uniformly 
in t G [0, T], This and (2.7) yield the desired result. □ 

2.3 Higher regularity 

Theorem 2.3. Let m be a nonneqative inteqer. We assume that q G D(A^ m+1 >' 2 ), h G 
D(A m / 2 ), and 

d'ff G C°([0, T]; D{A^ m ~ k ^ 2 )), k = 0,..., m. 

Then there exists a unique solution v to (2.1), such that 

d'fv G C°{[0,TY,D(A( m+1 - k V 2 )), fc = 0,l,...,m + l. (2.8) 

Moreover we have the estimate 

m +1 

E SUP ll 9 N(-> i )llwA(™+i-fe)/2) 

m \ 

E f\\c°([0,T]-,D(A(™~ k )/ 2 )) + llfl'llD(A( m+1 )/ 2 ) + IHI.D(A m /2) ] ■ ( 2 - 9 ) 

k =0 / 



(2.4) 

(2.5) 

( 2 . 6 ) 
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Proof, a) The proof is by an induction on m, the case m = 0 following from Proposition 2.1. 
b) We assume that the theorem is valid for some m £ N and suppose that 


f g £ D(A^ m+2 ^ 2 ), h £ R(H( m+1 )/ 2 ), 

\ a{ ! /€C ,o ([0,r];Z>(A(’ n + 1 - fc )/ 2 )), fc = 0,..., m + 1. 


We use the same strategy as in the proof of Proposotion 2.2. That is we differentiate (2.1) 
with respect to t and get that w := dtv is solution to (2.6). Next, using that h £ D(A^ m+1 ^ 2 ), 
/(0) -Age D(A m / 2 ) and 

d?(d t f) = d k+1 f e C°([0, T]; D(A^~ k V 2 )), k = 0,...,m, 

from (2.10), we get that d k+1 v = d k w e C°([0, T]; D{A ( e n+1 ~ k ' > / 2 )) for k = 0,1,..., m + 1, and 
the estimate: 

m+1 

sup ||9 t fc+1 u(-,t)|| jD(j4(m+ i-fc)/2 ) 
e[o ,T] 

m 

^ C ( E \\^t +1 f\\c°([0,T]-D(A<. m - k )/ 2 )) + ll^llo(A( m + 1 )/ 2 ) + 11/(0) — ^ffllz)(AW2) 

\ fc =0 




This entails £ C°([0, T]; H( J 4^ m+2 fc )/ 2 )) for k = 1,..., m + 2, and 

m+2 

E SU P ll 5 N(-^)ll£>UC-+2-fc)/2) 

^ *e[o.n 

/m+1 \ 

< c 1 E \\dtf\\c°([ 0 ,T]-,D(A( m + 1 - k )/ 2 )) + \\h\\D(A(™+ 1 )/ 2 ) + IMId(A(™+ 2 )/ 2 ) ] • ( 2 -H) 

\ fe =0 ) 

Further, as Av = f — d 2 v from the first line in (2.1), we find out that 

ll t ’(’> i )II.D(A( m +2)/2) < C (\\Av(-,t)\\D(A™/ 2 ) + IMI.D(A°)) 

< c (\\f('A)\\D(A™/ 2 ) + \\^t v {'A)\\D{A m / 2 ) + IMId(+>)) • (2-12) 

Here we used the identity ||u(-, t)\\ 2 D{A(m+2)/2) = II^M'^llElW 2 )+Efc=o \\A k / 2 v(-,t)\\ 2 D{A0) and 

the estimate P 1/2 u(-, t)\\ D{A0) < ELo P fe i>M)l \d(A°)- Since IMI^o) and \\d?v(-,t)\\ D (A™/ 2 ) 
are majorized by the right hand side of (2.11), uniformly in t € [0,T], (2.11)-(2.12) yield the 
assertion of the theorem for m+1. □ 


Remark 2.4. The result of Theorem 2.3 is similar to the one of [15, Sect. 7.2, Theorem 6], 
which holds for a bounded domain. This can be seen from the characterization of the D(A k / 2 ) 
for k = 0,..., m, displayed in Subsection 2.4. Namely, it is worth noticing that the m th -order 
compatibility conditions [15, Sect. 7.2, Eq. (62)] imposed on /, g and h, are actually hidden 
in the operatorial formulation of Theorem 2.3. 
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2.4 Characterizing the domain of A m l 2 for m E N* 

2.4.1 Elliptic boundary regularity 

In this subsection we extend the classical elliptic boundary regularity result for the operator 
V • cV , which is well known in any sufficiently smooth bounded subdomain of R n (see e.g. 
[15, Sect. 6.3, Theorem 5]), to the case of the infinite waveguide H under study. The proof 
of this result boils down to [24, Lemma 2.4] which claims elliptic boundary regularity for the 
Dirichlet Laplacian in f l. 


Lemma 2.5. Let r be a nonnegative integer. We assume that du is C r+2 and that c G 
W T+1 >°°(n) obeys (1.2). Then, for any ip G H r (Ll), there exists a unique solution v G H r+2 (Ll) 
to the boundary problem 


J —V • cVv = ip in Ll 
( v = 0 on dLl. 


(2.13) 


Moreover we have the estimate 

||^||r+2,r2 ^ C r ||^||r,Q? 


(2.14) 


where C r is a positive constant depending only on r, uj, the constant c m appearing in (1.2) and 

||c||jyr+l,°o(Q). 


Proof. The proof is by induction on r. 

a) We first consider the case r = 0. Due to (1.2) there is a unique solution v G Hq(Q) to (2.13) 
by the Lax-Milgram theorem. Moreover v satisfies the energy estimate 


w||i,n < C'||(/?||o ) f2, 


(2.15) 


where the constant C > 0 depends only on uj and c m . Here we used (1.2) and the Poincare 
inequality, which holds true in since uj is bounded. Furthermore, v is solution to the boundary 
value problem 


f — Av = f in Q 
\ v = 0 on <9fl, 


(2.16) 


where 

/ : = c _1 (^ + Vc- Vu). (2.17) 

Since / G L 2 (Q) then v G H 2 (Q) by [24, Lemma 2.4], and ||u|| 2 . q is upper bounded, up to 
some multiplicative constant depending only on oj, by ||/||o,n- As a consequence we have 


Mh.n < C' (Iloilo,a + ||«||i,n), 


from (2.17), the constant C' > 0 depending only on uj, c m and ||c|| w i.oo(q). This and (2.15) 
yield (2.14) with r = 0. 

b) Suppose that the statement of the lemma is true for r G N fixed, and assume that dco is 
C r+3 , c G W r+2 ’°° (Q) and ip G H r+1 (Q). Hence the solution v to (2.16) belongs to H r+2 {Tt) 
and satisfies the estimate (2.14), by induction assumption, and we have / G H r+1 (Q ) in virtue 
of (2.17). Further, v being solution to (2.16) where the boundary dco is C r+3 then v G H r+3 (Ll) 


8 


by [24, Lemma 2.4]. Moreover H^Hr+s^ is upper bounded (up to some multiplicative constant 
depending only on r and oj) by ||/|| r +i,fU From this and (2.17) then follows that 

IMIh-3,12 < C" (||<^||r+l,n + ||«||r+ 2 ,n) ) 

where the constant C" > 0 depends only on r, oj, c m and ||c||iyr+ 2 ,oo(Q). Putting this together 
with (2.14), we obtain (2.14) where r is replaced by r + 1, proving that the statement of the 
lemma remains valid upon substituting r + 1 for r. □ 

2.4.2 The domain of A m / 2 for m G N* 

In this subsection we characterize the domain of A m / 2 for rri £ N*. 

Proposition 2.6. Let m G N* and let k be either m — 1/2 or m. Assume that dev is C 2k and 
that c G P/ 2fc_1 ’ 00 (fl) fullfills (1.2). Then we have 

D(A k ) = {u G H 2k {Ll), u,Au, ..., A m_1 u G 7/ ( J(n)}. 

Moreover, the norm associated with D{A k ) is equivalent to the usual one in H 2k (Ll): we 
may find a constant c(k ) > 1, depending only on k, u, the constant c m defined in (1.2) and 
||c||w 2 fe-i,«>(s- 2 ), such that we have 

cikY 1 \\u\\ d{ a^ < \\ u hk,n < c{k)\\u\\ D{Ak) , u G D(A k ). 

Proof. It suffices to show that 

D(A k ) c {u G H 2k (Ll), u, Au, ..., A m ~ 1 u G Hq(Q)}, (2.18) 

and 

IM| 2 fc,n < c{,k)\\u\\ D ( A k), u G D{A k ). (2.19) 

The proof is by induction on m. 

a) We start with m = 1 and notice from the very definition of A 1 / 2 that D(A 1 ^ 2 ) = D(qy\) = 

Moreover we have 

P 1/2 «llo,n = Qa[u] > CmllVullo^, u G D(A 1/2 ), 

_i /o 

in virtue of (1.2). Bearing in mind that c m G (0,1), we obtain that ||u||i,n < c m 7 ||'w|lx?(^ 4 i/ 2 ) 
for every u G D{A 1 / 2 ). This establishes (2.19) for k = 1/2. 

Similarly, bearing in mind that D(A) = {u G Hq(Q), Au G L 2 (ll)}, we apply Lemma 2.5 
with r = 0 and ip = Au, where u G D(A) is arbitrary. We find that u G H 2 (Ll) satisfies 
IM| 2 ,n < Coll-Aullo^, which entails (2.18)-(2.19) for k = 1. 

b) Let us now suppose that the statement of the lemma is true for some m G N* fixed. Pick 
k G {m — 1/2, m} and assume that diw is C 2 ( fc+1 ) and that c G kF 2fc+1,00 (ll) satisfies (1.2). 
As D(A k+1 ) = {ri G D{A k ), Au G D(A k )}, we deduce from the induction assumption that we 
have 

D(A k+1 ) = {u G H 2k (Ll), Au G H 2k (Q) and u ,..., A m u G 
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with ||^4 J w||2fc,Q < cityWA 3 u\\ D (A k ) f° r 3 = 0,1. Thus, applying Lemma 2.5, with r = 2k and 
ip = Au, for u € D(A k+1 ), we get that u € H 2 ( k+1 \Q), proving (2.18) where (k + l,m + 1) is 
substituted for ( k,m ). Moreover, it holds true that 

IMl2(fc+i),n < C 2 k\\Au\\2k,n, 

and since Au E D(A k ), the induction assumption yields ||^4n||2fc,fi < c {k)\\Au\\ D (A k )- Therefore 
||n|| 2 (fc + i) j n is majorized, up to a multiplicative constant depending only on u, c and m, by 
IMIr> 04 fc +i), w hich is (2.19) where k + 1 is substituted for k. □ 

In view of Theorem 2.3 and Proposition 2.6 we obtain the following result. 

Corollary 2.7. Let m be a natural number. Assume that dui is C m+1 , that c € W m ’°°(n ) 
fulfills (1.2) and that (0 q,6i) E D(+ m+1 )/ 2 ) x D(A m / 2 ). Then the initial boundary value 
problem (1.1) admits a unique solution 

m+1 

tie p| C k ([0 ,T];^ m+1 - fc (H)). 

k= 0 


Moreover, we have 


m+1 

IMIc*([0,T];tf”*+l-fc(n)) < C (||0 O || m +l,n + |+|| m,£l) 5 (2.20) 

k=§ 

for some constant C > 0 depending only on T, ui and ||c||wm,oo/Q). 

3 Global Carleman estimate for hyperbolic equations in cylin¬ 
drical domains 

In this section we establish a global Carleman estimate for the system (1.1). To this purpose 
we start by time-symmetrizing the solution u of (1.1). Namely, we put 

u(x, t ) := u(x, — t ), i£ll, t G (—T, 0). (3-1) 

Under the conditions of Theorem 1.1, and since 0\ = 0, it is not hard to check that 

4 

tiG p| C k {[-T,T];H 5 ~ k (n)). 

k =3 

With a slight abuse of notations we put Q := Q x (—T, T), £ := T x (— T, T) and Qi := 
Pi x (—T, T), £l := <9ca x (—L, L) x (—T, T) for any L > 0, in the remaining part of this text. 
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3.1 The case of second order hyperbolic operators 

In view of establishing a Carleman estimate for the operator 

A:=A(x,t,d) = d 2 - V-cV + i?, (3.2) 

where R is a first-order partial differential operator with L°°(Q ) coefficients, we define for 
every <5 > 0 and 7 > 0 the following weight functions: 

if(x,t) = tpg(x,t) := \x' — <5a , | 2 — x 2 — f 2 and tp(x,t) = <pg,y(x,t ) := e 7 ^’^, (x,t) £ Q. (3.3) 

Further for L > 0 and T > 0 we introduce the space 

Xl,t := {m e H 2 (Q l ); ^ n u(-, ±L, ■) = 0 in w x (-T, T), d}u(-, ±T) = 0 in Q L , j = 0,1 j , 

(3.4) 

and state the: 

Proposition 3.1. Let A be defined by (3.2), where c obeys (1.2)-(1.7), and let l be positive. 
Then there exist do > 0 and 70 > 0 such that for all 6 > do and 7 > 70 , we may find L > £, 
T > 0 and sq > 0 for which the estimate 


s ]T s 2 ^\\e s ^ j Xjt v\\ 2 o,Q L <C 
3=0,1 


\\^ p Av\\l QL +s'£ t 
3=0,1 


. 2(1 ~j) 




lo ,8Ql 


(3.5) 


holds for any s > s 0 and v £ Xl,t- Here C is a positive constant depending only on u, a', ao, 
do, 7 o, so, Cm and cm- 

Moreover there exists a constant dp > 0, depending only on uj, t, do and 70 , such that the 
weight function ip defined by (3.3) satisfies 

<p(x',x n ,0) > d e , (x',x n ) £ux [-£,£], (3.6) 

and we may find e £ (0, {L — £)/ 2) and v > 0 so small that we have: 

max (p(x',x n ,t) < d£ := dee~ lu , |f| £ [T — 2e, T], (3-7) 

x£ujx [— L,L] 

max <p(x',x n ,t) < de, \x n \ £ [L — 2e, L\. (3-8) 

(x',t)£Z3x [—T,T] 


Proof. The proof is divided in three parts and essentially boils down to [20, theorem 3.2.1’]. 

First part: Definition of 6 o, L andT. Bearing in mind that — da, / 1 2 — Ij/' — da, / 1 2 = |^ , | — |y , | 2 — 
25a'-(x' — y') for all x',y' £ lo, we see that sup^/g^ \x'—da' 7 (x'—da '| 2 < |cu |(|o^|H-4d|oE, / 1) 

for every d > 0, where |w| := sup x ., ezJ \x'\. Hence the function 

f \ 1/2 

ge(6) = ( sup \x' — 5a | 2 — inf \x — 5a | 2 + £ 2 ) (3-9) 

\x'Guj x'Guj J 
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scales at most like d 1 / 2 , proving that there exists do > 0 so large that 

<5a 0 > ^^1 + —— Slifi) + M + 2 ^ cm + 2, d > <5o■ (3.10) 

Further, since ui is bounded and a' / 0 R n-i by (1.7), we may as well assume upon possibly 

1 /2 

enlarging do, that we have in addition Cm inf x / eZ j |x' — da'| > <+d) for all d > do- This and 
(3.10) yield that there exists i? > 0 so small that the two following inequalities 


dao — ( L + |w| + 2(1 + 


and 


x 'Go; 


Cm l/ 2 )) CM 2>0 ’ 

(3.11) 

da'| > T, 

(3.12) 


hold simultaneously for every L and T in (ge(S),gg(5) + d), uniformly in d > do- 


Second part: Proof of (3.5). We first introduce the following notations, we shall use in the 
remaining part of the proof. For notational simplicity we put x := (x,t) for (x,t) £ Ql and 
V x = (pi,... ,d n ,d n+1 ) T . We also write £' = ( 6 , • • ■ , U-i) T G ^ n_1 , £ = ( 6 , • • •, £n) T G M n 
and £ = (£i,..., £ n , £ n +i) T £ R n+1 . We call A 2 the principal part of the operator A, that is 
A 2 = A 2 (x, d) = df — c(x) A, and denote its symbol by A 2 (x,£) = c(x)|£| : — £ 2 +1 , where |£| = 

(E”=i • Since A 2 (x, V x +x)) = 4 (c(x)(|x' - da'| 2 + x 2 ) - x 2 n+1 ) for every x E Q L , we 

have 

+ (+V x +x)) > 0, x E Q l , (3.13) 

by (1.7) and (3.12). For all x E Q L and f £ M" +1 , put 


A 7 V' ^2 a a / 1 V' 

J(x,f) = J = ^ _ W + + 




j,k =1 


<9fc 




<90 / 




dfjdf k 


djf >, (3.14) 


where, for the sake of shortness, we write <+ j E N* +1 , instead of d/dxj, and x n+ i stands for 
t. Assuming that 

A 2 (x,i) = c(xm 2 ~e n+ i = 0, xeH, £ £ M n+ 1 \{ 0 }, (3.15) 


and 


V+ 2 (x, i) ■ V x +x) = 4 [c(x)(f' • (x' - da') - £ n x n ) + £ n+ ix n+ i] = 0, x E Q L , f £ M" + 1 \{0}, 

(3-16) 

we shall prove that J(x, £) > 0 for any (x, £) E Ql x M n+1 . To this end we notice that the 
first sum in the rhs of (3.14) reads (Hess(^>)V|.A 2 , V+ 2 ) = 8 (c 2 (|£'| 2 — £ 2 ) — £ 2 +1 ), and that 


71+1 

£ 

j,fc=i 


„ 3A 2 \ 0 A 2 <9 2 + 

4— ) -5^- - ( dkA )- 


dfj J df k 


dfjdf k 


°ji'= Y 


j,k=l 


ch4 2 \ dA 2 . . <9 2 A 2 
4— ) - 5 ^- - (d k A 2 )- 


dfj J df k 


dfM k 


djif 


2c Y ((++ _ i£i 

j,fc=i 


2 dfj 

d£k 


d k c) 9+ = 2 (2(Vc • £)(V^ • 0 - (Vc • V+|+) 
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since 


— (<9fc^2) ag 2 g| fc = 0 if either j or k is equal to n + 1. Therefore we have 

J = 2[4(c 2 (|^| 2 -a-d + i)+2c(Vc-0(V^-0-c(Vc-V^)|C| 2 ] 

= 4 [2c 2 (|£ , | 2 - £ 2 ) - (2 + (a/ - 5a') ■ V x >c - x n d n c ) £ 2 +1 - 2x n+1 £, n+1 Vc ■ f] , 

from (3.15)-(3.16). Further, in view of (3.15) we have c 2 (|£'| 2 — £ 2 ) > — c 2 |£| : > —c£ 2 +1 and 
|Vc ■ £| < |Vc||C| < (IVcI/c 1 / 2 )!^!, whence 


J > 4 


5a' ■ V x /c 


• Vc + 2c + 2 T 


|Vc| 

772 




(3.17) 


Here we used that fact that x n+ \ = t £ [0, T]. Due to (1.6)-(1.7), the rhs of (3.17) is lower 
bounded, up to the multiplicative constant 4£ 2 +1 , by the lhs of (3.11). Since f n +i is non zero 
by (1.6) and (3.15), then we obtain J(x,£) > 0 for all (x, £) £ Qi x R n+1 . With reference to 
(3.13), we may apply [20, Theorem 3.2.1’], getting two constants so = £ 0 ( 7 ) > 0 and C > 0 
such that (3.5) holds for any s > sq and v £ H 2 (Ql). 


Third part: Proof of (3.6)-(3.8). First, (3.6) follows readily from (3.3), with dt := e 7 ^ and 
/3g := inf x i Gu |x' — 5a'\ 2 — (f 2 . Next, for v £ (0,$) arbitrarily fixed, we put 

L = T = gi(5) + v. (3.18) 


Notice for further reference from (3.9), (3.12) and (3.18), that we have 

2> (3. 19) 

Cm \ C m J 

since c m € (0,1), by assumption. Similarly, as T 2 > g^{5) 2 +v 2 = sup x ; ew |a/ — 5a'\ 2 — ((3 2 — v 2 ), 
we dedude from (3.3) that 


<p(x',x n ,±T) < e y (. swp x'£“\ x '- Sa '\ 2 - x *- T2 ) < e^~ u2 )e-^, (x',x n ) £Ux [—L, L\. 

With reference to (3.19) we may thus chose e £ (0, (L — l)/ 2) so small that 

ip(x',x n ,t) < e~ lx2 , (x', x n ) £ ZJ x [-L, L\, \t\ £ [T - 2e,T], 

which entails (3.7). Finally, since t and x n play symmetric roles in (3.3), and since T = L, we 
obtain (3.8) by substituting (T, t) for ( L,x n ) in (3.7). □ 


3.2 A Carleman estimate for the system (1.1) 

In this subsection we derive from Proposition 3.1 a global Carleman estimate for the solution 
to the boundary value problem 


( d 2 u — V • cVrt = / in Q 
\ u = 0 on S, 


(3.20) 
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where / £ L 2 (Q). To this purpose we introduce a cut-off function x €= C 2 (M; [0,1]), such that 


xOn) : = 


1 if \x n \ < L — 2e, 


(3.21) 


0 if \x n \ > L — e, 

where e is the same as in Proposition 3.1, and we set 

u x (x,t) : = x(x)u(x,t) and f x {x,t) := x(x)f(x,t), (x,t) £ Q. 

Corollary 3.2. Let f £ L 2 (Q). Then, under the conditions of Proposition 3.1, there exist two 
constants s* > 0 and C > 0, depending only on u, t, Mq, t]q, a!, do, c m and cm, such that the 
estimate 


s £ s 2 ^>\\e^Vi,u\\l QL < C | lle^/H^ + S V 

3=0,1 


+ s S ~ 
3=0,1 




holds for any solution u to (3.20), uniformly in s > s*. 

Proof. Since u is solution to (3.20) we have 

d 2 u x — V • cVm x = f x + R\u in Q, 


where 

Ri = Ri(x,d) := [x, V • cV] = -(cAy + Vc • Vy + 2cVy • V), (3.22) 

is a first-order differential operator. Therefore, the function v(x,t) := rj(t)u x (x,t), where 
7] £ C 2 (M; [0,1]) is such that 


y(t) : = 


1 if \t\ < T — 2e, 
0 if |t| > T — e, 


satisfies 

d 2 v — V • cVv = g := ijf x + r/Riu + rfu x + 2?/3 t u x in Q. 

Moreover, as v(-,PL, •) = d Xn v{-, PL, ■) = 0 in u x (—T, T) and u(-,±T) = dtv(-,PT) = 0 in 
LIl, we may apply Proposition 3.1, getting 


» E S 2(1 - ) ' , ||e-«-vj it »|li, Ql < C 

3=0,1 



+ S 


E' 

3=0,1 


,2(1-3) 


\\^Kt<dQ L 


(3.23) 


Further, bearing in mind that 8Ql = 'Ll U (uj x {PL} x (—T, T)) U (Ql x {±T}), we deduce 
from the vanishing of v(-, PL, •) and V x ,tu( - , ±T, ■) in cu x (—T, T), and the one of u(-, ±T) and 
V Xi tv(-,±T) in LIl, that 


W^i,Ao,dQ L = lle^V^llo.E,, j = o, 1. 

Next we know from (3.8) and (3.22) that 

\\e s ' p r]Riu\\o i Q L < Ce kd< !MlL 2 (-T,T ; /P(n L ))> 


(3.24) 


(3.25) 
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and from (3.7) that 


\\e sip (ri"u x + 2r]'dtu x )\\o,Q L < Ce sde \\u x \\ H i^_ TT . L 2^ L ^. 
Hence, putting (3.23)-(3.26) together, we find out that 

« Y, » 2(1 -' , l|e*«’vj it o|| 1 Ql < C (l\e’ v f\\ 2 0 , QL + s *£ llo,^ + 

3=0,1 V 3=0,1 


(3.26) 


plsdt IL. 11 2 

e \\ u \\m(Q L ) 


(3.27) 


The next step of the proof involves noticing from (3.7) that \\G sip {l~r])'V :) xU x \\o^Q L < || V^7i- X ||o,Qz 

for j = 0,1, hence 


II^VXIkQ, < ||e^(l-r?)Vi % ||o, Qi + ||e^V^||o, Qi 

< e^||V> x || 0 ,Q L + ||e^V^|| 0 ,Q„ j = 0,1. 


(3.28) 


Furthermore, by combining the identity dtu x = (1 — rj)dtu x + rjdtu x = (1 — r))dtu x — rfdtu x +dtv 
with (3.7), we get that 

l|e S¥, 3 t w x || 0 ,Q L < \\e sip (l - v)d t u x \\o,Q L + ||e i <V x ||o,Q ll + ||e s ' p ^ t t;||o,o i 
< e sde (||<%u x || 0 ,Ql + 11 11 L°° (—T,T) ll^x llo,Qi,) + 
which, together with (3.28), yields 

E s2{1 ~ j) \\ eS ^it u x\\l,Q L <C E * 2(1_J ' ) + ||e^v>||^) . (3.29) 

j=o,i i=o,i 

Similarly, using (3.8), we derive from the identity d^u = d{u x + (1 — x)dt u f° r J = 0> 1 j that 
\\^d{u\\ 0:QL < \\e Sip diu x \\o,Q L + e sd t\\d J t u\\ 0 ,Q L , j = 0 , 1 , 
and from V x u = V x u x + (1 — x)V x ii — x / (0,..., 0, u) T , that 

\\e Sip V x u\\ 0) Q L < He^Vx^llo^i + e sde (||V x «|| 0j q l + \\x'\\l°°(-l,l)\H\o,q l ) ■ 

As a consequence we have 

E s^lle^Vi^llo,^ < C E s 2(W) 1,q l + \\e av Vi, t Ux\\o,Q L ) • (3-30) 

3 = 0,1 7 = 0,1 


Finally we obtain the desired result by gathering (3.27) and (3.29)-(3.30). 


□ 


4 Inverse problem 

4.1 Linearized inverse problem and preliminary estimate 

In this subsection we introduce the linearized inverse problem associated with (1.1) and relate 
the first Sobolev norm of the conductivity to some suitable initial condition of this boundary 
problem. 
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Namely, given c* £ Ap for i = 1,2, we note Ui the solution to (1.1) where q is substituted 
for c, suitably extended to (—T, 0) in accordance with (3.1). Thus, putting 

c := ci — C 2 and f c := V • (cVu 2 )j (4-1) 


it is clear from (1.1) that the function u := u± — U 2 is solution to the linearized system 

( df u - V • (ciVu) = f c in Q 

< u = 0 on £ (4.2) 

\ u(-, 0) = dtu(-, 0) = 0 in 12. 


By differentiating fc-times (4.2) with respect to t , for k £ N* fixed, we see that u^ 
solution to 


dfuW - V • (ciVnW) = / c (fc) in Q 
u^ = 0 on £, 


d£u is 
(4.3) 


with fc''- 1 := df f c = V • (cVu^), where u 2 k> stands for d^U 2 - 

We stick with the notations of Corollary 3.2. In particular, for any function v, we denote 
XV by v x , where x is defined in (3.21). Upon multiplying both sides of the identity (4.3) by 
X, we obtain that 


dt u x k) ~ v ' ( c i v «x fc) ) = /£> - 9k hr Q 
u x ^ =0 on E, 


(4.4) 


with 


fc x ■= V • (c x Vn 2 ) and g k := V • (ci(V X )u (fe) ) + ciV X ■ W fc) + cV X ■ Vu^\ (4.5) 

Notice that g k is supported in 1 i e := {x = (x',x n ), x' £ w and \x n \ £ (L — 2e,L — e)}. 

Having said that we may now upper bound, up to suitable additive and multiplicative 
constants, the e S¥, i' ,0 i-weighted first Sobolev norm of the conductivity c x in 12^, by the corre- 
sponding norm of the initial condition u x (-,0). 

Lemma 4.1. Let u be the solution to the linearized problem (4.2) and let x be defined by 
(3.21). Then there exist two constants s* > 0 and C > 0, depending only on uj, e and the 
constant Mq defined by (1.12), such that the estimate 

£ < Co' 1 ( £ || e «o<-.»>v>»W(-,0)g, ni + e 2«fc 

l=o,i \i=o,i 

holds for all s > s*. 

Proof. Let 12* be an open subset of with C 2 boundary, such that 

w* x (—(L — e), L — e) C $2* C w* x (—L, L), (4.6) 

where e is defined by Proposition 3.1. We notice from (3.21) and (1.11) that djc x € 1/q( 12*) 
for all i £ N* and j = 0,1. 
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Further, with reference to (1.10) we may assume upon possibly enlarging 5 £ [doj+oo); 
where do is the same as in Proposition 3.1, that we have 

|V0 O • (xi - Sai,... ,x n -i - 5a n - i, -x n )\ > > 0, x £ 11*. 

Thus applying [18, Proposition 2.2] 5 with D = SI*, P(x,d)v = V • ((y0o)v) and v = djc x £ 
Hq(Q*) since x( x n) = 0 for x n > L — e, for i £ N* and j = 0,1, we obtain that 

S || e ^’°)^c x || 2 A < C||e^-’°)V • ((dic x )Ve 0 )||g,n., * G N*. j = 0,1. (4.7) 

Since c x (x',x n ) = 0 a.e. in w* x ((— L,— (L — e)) U (L — e,L)) by (3.21), we have 
||e S¥, (' ,0 )^c x ||o,f 2 , = ||e S¥ ’^ ,0 ^"c x ||o,n I/ for each i £ N* and j = 0,1, from (4.6). We derive from 
this and (4.7) that 

a ||e^ ( -- 0) ^c x ||g ini < C||e^’°)V • (O/c x )W 0 )|liW, i € N*, j = 0,1. 

Further, taking t = 0 in the first line of (4.4) with j = 0, we get that 

V • (c x W 0 ) = 4 2) (‘> °) + cV X ' V0 O . 

From this, (4.8) and (3.8) then follows that 

«l|e s ' p( -’ 0) c x ||g,n i < c(||e^-’ o )4 2 )(-,O)||^ i + || e ^- o ) c Vy.V0o||^ i 

< c-(|| e s ^°)4 2 )(-,o)|| 2 >ni + e 2 ^). 

Similarly, since V • ((<9jC x )V#o) = <9iV ■ (c x 0 o) — V • (c x V<9j0o) for every i £ N*, we derive from 
(4.9) that 

v • ({d iCx )ve 0 ) = d iU W{; o) + di{cve 0 • v x ) - Vc x ■ vd t e 0 - c x Ac^ 0 . 

As a consequence we have, 

|| e ^(-,°)V. ((a i c x )V0 o )||^ 


(4.8) 

(4.9) 


(4.10) 


< Cl ^||e^-’°)^4 2 )(-,0)|| 2 A + ||e^’°Vc x || 2 A + e 2 ^J , i G N*, 

according to (3.8). Summing up the above estimate over i in N*, it follows from (4.8) that 


s\\e 


s<f(-,0) 


Vc, 


2 

lo,fi L 


< C ( ||e^(-.°)v4 2 )(-,0)|| 2 nL + ]T ll^ ( '’ 0) V J c x || 2 Oi +e 

l=o,i 


2 sdp 


This and (4.10) yield the desired result. □ 

s Let D be a bounded open subset of R n , n > 1, with C 2 boundary, and consider the first-order operator 
P{x,d) := J2i=iPi( x )di +Po(x), where po € C°(D) and p := (pi,... ,p n ) € C 1 (D) n . Assume that 

| p(x) ■ (an - 5a i, ... ,x n -i - 5a„-i,-x n )\ > p m , x £ D. 

for some p m > 0. Then for any p m > max{||po|lc«(l 3 )i llPdlcifo)" * e N*}, there exist s* > 0 and C > 0, 
depending only on pM, such that the estimate 

l|e , ' p( ' ,0) «llo,D <Ca- 1 ||e^-' 0) Pt;||g iD , 

holds for all s > s, and v € Hq(D). 
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4.2 Completion of the proof 

The proof is divided into three steps. 

Step 1. The first step of the proof is to upper bound u x \-, 0) in the e S¥, ^’°hweighted H 1 (£Il)- 

(2) (3) 

norm topology, by the corresponding norms of u x J and u x ; in Ql, with the aid of the following 
technical result, borrowed from [3] [Lemma 3.2]. 

Lemma 4.2. There exists a constant s* > 0 depending only on T such that we have 

I \ z i'i 0)llo,n L ^ 2 ( s ll z llo,Qi + s l \\ d tA\l,Q L ) j 
for all s > s* and z G H l {—T, T; L 2 {TIl))- 

Namely, we apply Lemma 4.2 with z = e stp dju x ' for i € N* and j = 0,1, getting 
|| e M-.0) fl |„P) ( . i O)||2 ni < 2 ( S || e >^»«||5 i(3i +S - I || e .^„(3)||2 gi ) , „ > 

Summing up the above estimate over i G N* and j = 0,1, we obtain for all s > s* that 

£ ||e^V£ 2 )(-,0)||£ L < 2 £ ^s\\e s ^u^\\l QL +s-^e s ^u^^\l QL ) . (4.11) 

l=o,i l=o,i 


Step 2. The next step involves majorizing the right hand side of (4.11) with 

Ms) := £ ^MUe^Vi^llo.SL. k = 2 > 3 - 

i=o,i 


(4.12) 


(k) 

Indeed, since , for k = 2,3, is solution to (3.20) with c = c\ and / = f c — according to 
(4.4), then Corollary 3.2 yields 

S E S 2(1 - j) ||e*»vj il aW||S, ei < C (||e^/ij>||g, Qi + ||e**> 9 *||V + S 3 e 2 **||«W||?, gi + 


1 = 0,1 


for s large enough. In light of (4.11) this entails that 

E ||e‘*(' o) V%«(-,0)||§, ni 


1 = 0,1 


< 


c E (l|e w /S > llo, 0 I . +11^112,4+« 3eM ‘ll“? ) lli,«i.+«(>*(« 

k= 2,3 


(4.13) 


Further, recalling (4.5), we see from (1.5) and (1.12) (resp., from (1.5), (1.6), (1.12) and (3.8)) 
that the first (resp., second) term of the sum in the right hand side of (4.13) is upper bounded 
up to some multiplicative constant, by £ =0 i ||o S¥: ’V- 7 c x ||q q l (resp., e 2s ^(||u( fc ) || 2 + !))■ 

From this and Lemma 4.1 then follows for s sufficiently large that 


0 E l|e“'’ l ' , " , VEllo,n i 
1 = 0,1 


< C| E \\e“ p V 1 c x \\l QL +e 2 ‘ d ‘+ E (»V to *||«“>||l gi +»W» 

0=0,1 k= 2,3 


• (4-14) 
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(4.15) 


Step 3. Finally, we notice from (3.3) that 

l|e* y V>c x ||„, Ql = i = 0,1, 

where p s (x) := fl e i ' i: e ,: '' : ' f .' 'reedii f [t = J 2 ), e rwf-rfw' cp: ,t 2 ))^ f or a ]] a ^ Bearing in 
mind that <^(x,0) > n := e 7 ^ mf ^'6“l x,_5a, l _Z/2 ^ > 0 for all x € we get that 

\\Ps\\l^ l ) < j\- 2s ^- e ^-^dt, s > 0. 

Therefore we have lim^+oo p s = 0, uniformly in by the dominated convergence theorem, 
so we derive from (4.14)-(4.15) that 

a V || e M-.0) V J Cx g ni < C ( e 2 ** + £ ( s V**||B<‘)||?, gi + shr(a)) ) . (4.16) 

3 = 0.1 \ k= 2,3 / 

upon taking s sufficiently large. With reference to (3.6)—(3.8), this entails that 

E llVElliW < C E (s 2 e- 2s ^~^ + b *(*)) . (4.17) 

. 7 = 0,1 k= 2,3 

Here we used ( 1.5)-( 1.12) and the embedding Qg C f2^ in order to substitute for in the 
left hand side of (4.16). Now, taking into account that dp < dg, we end up getting the desired 
result from (4.17). 
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